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The subject of this paper is to present a theory for the dynamic radiation force produced by 
dual-frequency ultrasound beams in lossless and non-dispersive fluids. A formula for the dynamic 
radiation force exerted on a three-dimensional object by a dual-frequency beam is obtained stemming 
from the fluid dynamics equations. Dependence of the dynamic radiation force to nonlinear effects 
of the medium is analyzed. The developed theory is applied to calculate the dynamic radiation force 
exerted on solid elastic (brass and stainless steel) spheres by a low-amplitude dual-frequency plane 
wave. The dynamic radiation force is calculated by solving the dual-frequency plane wave scattering 
problem for solid elastic spheres. Results show that the dynamic radiation force on the analized 
spheres presents fluctuations similar to those present in the static radiation force. Furthermore, 
analysis of the static and the dynamic radiation forces on the stainless steel sphere shows that they 
have approximately the same magnitude. 

PACS numbers: 43.25.+y, 43.35.+d 



I. INTRODUCTION 

It is known that acoustic waves carry momentum. When an acoustic wave strikes an object jpart of its momentum 
is transferred to the object giving rise to the acoustic radiation force phenomenon 0,00,01. Acoustic radiation 
force has found practical importance in many applications. For example, measure thepower output of transducers 
in medical ultrasound machines Q, ultrasound radiometer @|, liquid drops levitation Q, and motion of gas bubbles 
in hquids Q. 

In these apphcations, the radiation force is static, being generated by a continuous-wave (CW) ultrasound beam. 
Time-dependent (dynamic) ultrasound radiation force can also be produced by amplitude-modulated (AM) or pulsed 
ultrasound beams. In general, an AM beam produces a harmonic radiation force at the modulation frequency, while 
a pulsed beam generates a pulsed radiation force. Dynamic radiation force has been applied for measuring the 
ultrasound power of transducers using a disk target Q or a shaped- wedge vane ) and determining ultrasound 
absorption coefficient in liquids [Tl]|. 

In recent years, dynamic ultrasound radiation force has become of practical importance in elastography; specifi- 
cahy in the fohowing imaging techniques: 



Acoustic radiation force impulse imaging (ARFI) uses pulsed ultrasound radiation force to produce displace- 
ment in tissue which is detected to form an image of the tissue • 

Shear wave elasticity imaging (SWEI) images tissue properties by detecting shear acoustic waves induced by 
the harmonic ultrasound radiation force produced by an AM ultrasound beam 

Vibro-acoustography maps the mechanical response of an object to a harmonic ultrasound radiation force 
produced by two overlapping CW ultrasound beams whose frequencies are slightly different In this 

context, it has been demonstrated that viscoelastic properties of gel phantoms can be accurately determined 
by measuring the amplitude of vibration induced by the dynamic ultrasound radiation force on a small sphere 
embedded in the medium llfll. 
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FIG. 1: Diagram showing the theory of acoustic radiation force in different propagating mediums. Dashed ellipses show the 
lack of theory for acoustic radiation force. 

Lord Rayleigh was the first to propose a theory for acoustic radiation force in lossless fiuids due to compressional 
waves Static radiation force in lossy fluids has been studied by Jiang et al. Doinikov ^19|, and Danilov 
et al. j2y| A study of the static radiation force in lossless isotropic elastic solid can be seen in Ref [2jj. Usually 
the radiation force exerted on an object target by a CW ultrasound can be obtained by solving the vector surface 
integral of the radiation-stress tensor (defined as the time-average of the wave momentum fiux) over a surface 
enclosing the object. The radiation-stress is obtained in terms of the incident beam and the scattered field by the 
object. Several authors derived the static radiation force, by solving the scattering problem of CW plane waves by 
spherical or cylindrical objects (23. 123. 12R . Most studies realized for the dynamic radiation force 

have been focused on finding applications to it. No theoretical efforts to tackle the problem of dynamic ultrasound 
radiation force exerted on an embedded object in a medium have been made whatsoever. Figure depicts the 
theoretical realm of acoustic radiation force. This figure includes the contribution of this paper: dynamic radiation 
force in lossless fiuids. Dashed ellipses show the lack of theoretical models for acoustic radiation force. 

The increasing applications of dynamic radiation force of ultrasound in elastography techniques prompted us to 
develop a method to calculate this force. Here, we present a theory of dynamic ultrasound radiation force exerted 
on arbitrary shaped three-dimensional objects. The theory is restricted to the radiation force produce by dual- 
frequency CW ultrasound waves in lossless fiuids. In what follows, we briefiy discuss the dynamic equations of 
lossless fiuids in Sec. Ill Al In Sec. IIIBI we present a theory of ultrasound radiation force. We obtain a formula 
for the dynamic ultrasound radiation force exerted on an object by a dual-frequency CW ultrasound wave. The 
formula is given in terms of a vector surface integral of the wave velocity potential over the object's surface. Explicit 
dependence of the dynamic radiation force with the medium nonlinearity is pointed out. In Sec. IIIII we apply the 
theory to calculate the dynamic ultrasound radiation force exerted on solid elastic (brass and stainless steel) spheres. 
Finally, in Sec.lTvl we summarize the main results of this paper. 

II. THEORY 

A. Dynamics of lossless fluids 

Consider a homogeneous isotropic fiuid in which thermal conductivity and viscosity are neglected. This corre- 
sponds to the so-called ideal fiuid. The medium is characterized by the following acoustic fields: pressure p, density 
p, and particle velocity v. Here all acoustic fields are functions of the position vector r G R'^ and time i S E. In 
an initial state without sound perturbation these quantities assume constant values given by p = po, p = po, and 
V = 0. The quantity p — po is the excess of pressure in the medium. The equations that describe the dynamic 
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of ideal fluids can be derived from conservation principles for mass, momentum, and thermodynamic equilibrium. 
These equations, neglecting effects of gravity, are presented here as follows (29i |: 

|^ + V-(pv)=0, (1) 

p'i - -v., (2) 

where V is the gradient operator, the symbol • stands for the scalar product, and B/A is the Fox-Wallace 
parameter [^^l which characterizes the nonlinearity of the fluid. Equations (C]|-(0I form a system of nonlinear 
partial differential equations that gives a full description of the wave propagation in the fluid. The conservation 
of fluid mass is represented by Eq. Q. Equation l(2l is a version of the Newton's second-law in fluid dynamics. 
Equation ijSl is an adiabatic equation of state known as Tait's equation. 

A lossless fluid is irrotational. Thus, according to the Helmholtz vector theorem, the particle velocity can be 
expressed in terms of the velocity potential function cf) as 

V = - V0. (4) 

The velocity potential can be expanded as a sum of a linear term and higher-order contributions as follows 

= + (5) 

where 0^^^ and 0^^^ are the linear and the second-order velocity potentials, respectively. In terms of the linear 
velocity potential, the linear pressure and velocity fields are given by 

v(i) = -V(/.(i). (7) 



B. Instantaneous net force 



A volume element in a fiuid is subject to a stress caused by the sound wave propagation throughout it. Stresses 
caused by sound perturbation in the fluid should be described by Eq. (j^J. 

Consider an ultrasound beam striking a homogeneous object of finite extension and surface Sq at rest. As the 
ultrasound field hits the object, its surface may be deformed and dislocated. We denote the object's surface at the 
time t by S. Fig. |2l depicts the interaction between the incident ultrasound wave and the object target. 

The instantaneous net force f acting on the object is obtained by integrating Eq. ^ on the object's volume. 
Since the ambient pressure po does not contribute to the net force on the object, we can replace p by p — po in 
Eq. Hence, integrating Eq. (EJ on the object's volume and using the Gauss' integral theorem, we obtain 

{ ^ - f {p - Po)ndS, (8) 
Js 

where n is the outward normal unit-vector of the integration surface. 

Radiation force is a phenomenon that depends on the interaction of second-order acoustic fields with the object 
target. To avoid the integration of p — po over the time-dependent object's surface S, we need to find p — po up to 
second-order. The idea is to replace S by Sq for second-order integrands in Eq. For first-order integrands, we 
should keep S and analyze the contribution of the integral to the radiation force. From Eqs. Q - (EJ, one can show 
that the second-order excess of pressure is given by 
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FIG. 2: Net force exerted on an object by an ultrasound beam. The dotted contour depicts changes on the object's surface. 



where Cq is the small-amplitude speed of sound. Substituting Eq. ^ into Eq. ||SJ and holding terms up to second- 
order, we find 



By using the relation 



Js dt dtJs"^ is 



in Eq. IjlOII and again keeping up to the second-order terms, we get 

d 



-Po- 



>^ S J So / So 



\-IS J So 

We call the quantity T as the radiation-stress tensor. It is given by 



n • JdS. 



(11) 



P 



(1)2 



Po(vW-vW) 



l + Pov«v«, 



(12) 



where I is the 3 x 3-unit matrix and pov'^W^^ is a dyad [131 known as the Reynolds' stress. We have written down 
Eq. Ijl2ll using the identities n • v'^^V^^^ = (n ■ v'^^)v^^-' and n • I = n. 

Consider g as a function of time. The Fourier transform of g is defined as 

r + oo 

g{co) = T[g] ^ / g{t)e-^-'dt, 



where u (angular frequency) is the reciprocal variable of t and j is the imaginary-unit. To analyze the frequency 
components present in the net of force, we take the Fourier transform of Eq. lfTT|l as follows 

J'[T]dS. (13) 

This equation gives a description of the frequency spectrum of the net of force acting on the object. We shall use 
this equation to calculate the static and the dynamic radiation forces. 
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C. static radiation force 



The static component of the net force f is commonly known as acoustic radiation force. Here, we call it static 
radiation force. Acoustic radiation force has been studied as a time averaged force. The rule of the time-average is 
to isolate the static component (w — 0) of the net force acting on the object. 
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Consider an incident ultrasound wave, periodic in time, striking the object target. The static component of the 
radiation force corresponds to w = in Eq. The first two integrals in the right-hand side of this equation 

become zero. Therefore, the static radiation force is 

f, = - / n • T[T]^^odS. 

J So 

Recognizing that the time-average of T over a long time interval is (T) = J'[T]ui=o, we get 

fs = - / n • {J)dS. (14) 

■J So 

Note that fs is a real quantity. 

The static radiation force can be understood as follows. The incident ultrasound beam is scattered by the object. 
The static radiation force is the time averaged rate of the momentum change due to the scattering by the object. 

The time-average of the radiation-stress is a zero-divergence quantity, i.e., V • (T) = 0|. This means that no 
static radiation force is present in an ideal fluid if there is no target. Consequently, steady streaming does not 
happen in lossless fluids |3^ . 



D. Dynamic ultrasound radiation force 

In this section, we concentrate on the dynamic radiation force produced by an amplitude-modulated (AM) ultra- 
sound beam. We chose the dynamic ultrasound force produced by AM ultrasound waves because of its importance 
in elastography; speciflcally in SWEI and vibro-acoustography. An AM ultrasound wave is described by a carrier 
/o and modulation A/ frequencies. Usually the carrier frequency is much larger than the modulation frequency. 

Here, the AM ultrasound beam is produced by two overlapping CW ultrasound beams. We call this beam a 
"dual-frequency ultrasound beam". The frequencies of this beam are /a = /o + A//2 and fb = fo — A//2. The 
corresponding angular frequencies of the beam are uja — 27r/a and ui, = 27r/fc. The frequencies /o and A/ are also 
called the center and the difference frequencies, respectively. 

When an object is placed in the wave path, the dual-frequency ultrasound beam is scattered by the object. The 
flrst-order velocity potential in the medium is given by 

^^^^ =^{4>ae'^^' + 4>be'^''}, (15) 

where is the real-part of a complex quantity. The functions (j)a and (j>b are spatial complex amplitudes of 

each frequency component of the wave. The amplitude functions in Eq. Ijl5|l are given as the sum of the velocity 
potential amplitudes of the incident and scattered waves. 

The radiation force exerted on the object by the dual-frequency beam includes the following components: a static 
component (w = 0), a component at the difference frequency Aw, and high-frequency components at 2(jJa, 2ujb, and 
oja + ujb. In this paper, the component at Auj is called the dynamic radiation force. 

To obtain the static and the dynamic radiation force produced by the dual-frequency ultrasound beam, we 
calculate the Fourier transform of the radiation stress T up to the difference frequency Auj. Substituting Eq. ((T5ll 
into Eq. lfT2|l . through Eqs. Q and and taking the Fourier transform of the result, we obtain the static, (T), 
and the dynamic, Taw, radiation-stresses as follows 

(T) = (T„) + (T,), (16) 

tA.. = y [[kah^ak - V4>a ' V^^) I + V^a^^ + ^ ^a] ■ (17) 

The quantities (Tq) and (Tfc) are the averaged radiation-stress of each frequency component of the dual-frequency 
wave and ka = uJa/co and kb — uob/co. ExpHcitly, the averaged radiation-stresses (Ta) and (Tfc) are given by 

(T„} = ^(fc^|<^™p-|V^™p)l + ^(V0„VC), m = a,b. (18) 

Now both the static and dynamic radiation force can be obtained. According to Eq. (|14|l the static radiation force 
is 

n • ((T,) + (Tb))d5. (19) 

So 
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Based on Eq. lfT3|l the amplitude (in complex notation) of the dynamic radiation force at Aw is 



-jpoAuj < T 



1(2) 



ndS 



So 



(20) 



So 



where 41 = ^[<I>^^'^]^=A^- 

Let us analyze the contribution of the first integral in the right-hand side of Eq. l(2?Hl . By using the Gauss' integral 
theorem, this integral can be expressed as 



Vo 



sv 



where Vq is volume of the object at rest and SV is the volume variation of the object induced by the excess of pressure 
at the time t (see Fig. |2ll . The first integral in this equation does not have any harmonic term at the difference 
frequency Aw, thus, by taking the Fourier transform of it and isolating the frequency component uj = Aw, we obtain 





[ 0(1) nd5 


= -T 


[ v(i)dy 






Js 


UJ—/S.UJ 


Jsv 


uj—A(jj 



Inasmuch as the particle velocity is a limited and continuous function inside the volume SV, there exists a point 



ro £ 6V such that | 
term in Eq. Ij2fl|l is 



JgyV^-'-'^dV = 5Vv'^^\ro,t). Thus, the amplitude of the dynamic radiation force of the first 



b'-^'^ndS 



^Au 



5M^r, 



d=Au 



(21) 



where Vq^' — v(i)(ro,t) and 5M — po5V corresponds to the fluid mass variation at the time t caused by the object 
vibration. Therefore, the dynamic radiation force can be written as 



_ Po 
2 



5Mv, 



(1) 



j=Au 



I So 



ndS 



{kah4>ak - V0a ' V0,*) n + (n • V0a)V0^ + (n • V^^V^a 



dS. 



(22) 



Thus, the dynamic radiation force exerted on the object target comes from three different interactions between 
the ultrasound fleld and the object. The flrst term in the right-hand side of Eq. ll22ll corresponds to momentum 
rate exchange due to fluid mass variation caused by the object vibration. The next term in this equation accounts 
for the interaction of the second-order velocity potential with the object target. None of these terms are present on 
the static radiation force formula itToll . The last term in Eq. (l22ll is related to the interaction of the radiation-stress 
and the object. We shall apply Eq. lf22|l to calculate the dynamic radiation force on a solid small sphere. 

In time-domain, the dynamic radiation force is given as the inverse Fourier transform of {auj- Therefore, the 
dynamic radiation force exerted on a object is given by 



fA^(t) = K{fA^e^^'"*} 



(23) 



Note that if Aw = and (j)a = (f>b, we have — is- Consequently, the magnitude of both dynamic and static 
radiation forces becomes equal when Aw = 0. 



III. DYNAMIC RADIATION FORCE ON A SOLID SPHERE 



A. Linear ultrasound scattering 



Consider a collimated dual-frequency plane wave described by Eq. itTfill impinging on a soHd elastic sphere of 
radius ro localized at the origin of the coordinate system. The beam propagates along the z-axis. The sphere 
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is characterized by three parameters: density pi, compressional wave speed Cc, and shear wave speed Cj. The 
amphtude functions of the incident plus the scattered waves are given, in spherical coordinates, by 



= Ao 5^(2n + l)(-jr b„(fc„r) + 5„,„/i(f^(fc„r)]P„(cos0), 



— a, b. 



(24) 



where Aq is the magnitude of the incident wave, km = ^m/co, jn is the spherical Bessel function, Sm,n is the 

(2) 

scattering function determined by boundary conditions, and hn is the spherical Hankel function of second-kind. 
The scattering function is given by 



^m^n^n (^m) ^m^n (^m) 



(25) 



where the prime symbol means derivative with respect to the function's argument and Xm = kmro- The coefRcient 
Dm,n is given by 



■Drn . n. 



2pi 



2n{n + l)jn{Xyn.s) 



(a^m.s/S - ^2 + n)j„(a:™_c) - '2x^^cjn+i{Xm,c) 2n{n + 1)[(1 - n)j„(x„^,j) + 2:™^si«+i(a:™,5)] 



(?2 -l)jn ('^m,c) ^rn.cjn-^-l i^Xm.c) 



-1 (26) 



where Xm,c = (co/cc)a;m and Xm,s = (co/cs)^^. 

The solution for the rigid and movable sphere can be obtained by setting Cc,Cs ^ oo in the previous discussion. 
The solution for liquid spheres is achieved by letting Cs — > 0. 

B. Second-order ultrasound scattering 

Before calculating the dynamic radiation force, we need to analyze the contribution of the second-order velocity 
potential in Eq. II22|I . The amplitude function for ^'^ incident wave is calculated in the Appendix. In the 
pre-shock wave range, it is given by 



/;(2) - 



ev. 



a_--]Akz 



2Auj 



(27) 



where e = l + B/ {2A), vq is the peak amplitude of the velocity potential at the ultrasound source, and Afc — Auj/cq. 
To simplify our analysis it was assumed that Akz ^ 1. If the difference frequency is about 10 kHz, then the target 
should be around 1 cm of the ultrasound source. 

Now, we need to solve the scattering problem for the second-order velocity potential. This problem is similar 
to the linear scattering presented in Sec. IIII Al Hence, the total second-order velocity potential amplitude can be 
written in spherical polar coordinates as 



(28) 



where the scattering function Sn is given through Eqs. Ij25ll and 112611 by setting Xm = Akrg. 

According to Eq. Ij22ll the contribution of the second-order velocity potential is given by integrating Eq. Ij28ll over 
the sphere surface. This contribution to the dynamic radiation force has only one component in the z-direction 
given by 



Fi = ~j2nrlpoALu / (^^'^(ro) sin cos = -nrlEo^MAkro) + ^ (Afcro)], 



(29) 



where Eq = {poVq)/2 is ultrasound energy at the wave source. 
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C. Dynamic radiation force function 



To calculate the dynamic radiation force we introduce the following variables Um — k^r (m — a, b) and v — coaO. 
By symmetry considerations the dynamic radiation force on the sphere is only in the z-direction. By substituting 
Eq. lf2ill into Eq. ll22ll leads to the amplitude of the dynamic radiation force as 



where the amplitude functions are 

F2 



F. = 



Fa 



-nrlpokah 
-irrlpokakb 



/Ac. = Fo + F^+F2 + F3 + Fa + F5 



1 dUa 



dub 



1 dv 



dv 



2\ 



F. = 



-Trropo 



vdv, 
v{l — v'^)dv, 
kh- 



dv 



dub 



(1 - v^)dv, 



(30) 

(31) 
(32) 
(33) 
(34) 



and Fq and Fi are given by Eqs. II21II and Ij29|l . respectively. 

To obtain the amplitude functions of the dynamic radiation force in Eq. (13011 . we substitute Eq. Ij24ll into Eqs. (I31|l - 
(PHl, which leads to 



Fo =- 



Fi 
F5 



Xa^b 



^{n + l)xaXb{R a, 71 



n+l 



I- Ra,n+lRb,n)' 

F3 = (" + '^)XaXbiRa,nR'b,n+l + R'a.n+lR'*b,n) ^ 

nr _ -7*7. ' ^ ' ' ' ' 

00 

^ n{n + \){n + 2){Ra,nRl„+i + i?a,«+ii?fc,„), 



ri=0 



XaXb 
XaXb 



71=0 



^[7l(n + l){xbRa.nR'ln+l + XaR'a,n+lRb,n) 



n=0 

n+l){7l + 2){XaK,nK 
,(2) 



(35) 
(36) 
(37) 

(38) 



n+l 



XbRa,n+lR' 



b,n)\j 



where i?„i,ri — {~j)"'[jn{xm) + Sm.nhn (xm)] and Eauj — PokakbAQ is the difference frequency component of the 
ultrasound energy density. 

Let us focus on the contribution of Fq to the dynamic radition force on the sphere. The magnitude of the velocity 
particle of the incident plane wave is po / (poCo) , where po is the magnitude of the incident pressure. The magnitude of 
the velocity particle inside the object volume variation SV is smaller than its counter-part in the fluid. Hence, from 
Eq. lf2T]l we have |Fo| < Aw(5M„iaxPo/(poCo), where SM^^^ is the maximum amount of fluid mass dislocated by the 
sphere vibration. Measurements of the amplitude dislocation induced by the dynamic radiation force (A/ < 1 kHz, 
/o around 1 MHz, and po < 60 kPa) on a stainless steel sphere of 1 mm diameter in water, yielded results less than 
1 pm [^j. Now, we compare the magnitude of Fq and F2. From Eq. (j35ll we have 

\Fo\ ^ Pof^AfrlSr 

\F2\ CoPo 

where 5r is the sphere radius variation. For the given condition of the experiment on measuring the dislocation 
amplitude of the sphere caused by the dynamic radiation force, the ratio |Fo|/|i^i| ^ 10~^. In this case or whenever 
Sr is negligible, only the components Fi to F5 are relevant to the radiation force formula ij^fijl . 
One can show from Eq. (I23|l that the dynamic radiation force on the sphere in time-domain is 



(39) 
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TABLE I: Physical constants used in the calculation of the radiation force functions. 





Speed of sound 


Material 


Density Compressional Shear 




Kg/m^ m/s m/s 


Brass 


8100 3830 2050 


Stainless steel 


7670 5240 2978 



where Y^^i is the dynamic radiation force function given by 

2 °° 

Yauj = X^^"" ^)[iXaXb - n{n + 2)){Ra,nRl„+l + Ra.n+lRln) 

^"■^'^ n=0 

+ n{xhRa,nR' b,n+l + ^aR'a.n+lRljn) ~ {n + 2){XaR'a.nRb,n+l + XbRa,n+lR' b,n) 

+ XaXb{R'a^nR'l,n + l + R'a,n+lR' b.n)] - Rq, (40) 

where Ro = ||f^[ji(Afcro) + Sih^^\Akro)]. 

The sphere is also subjected to a static radiation force, which is the sum of the force due to each ultrasound wave 
in the incident beam. The static radiation force on a spherical target has been calculated by Hasegawa j23|. The 
result reads 

fs=Trrl{EaYa + EbYb)e,, 
where Em — po(fcm^o)^/2 (rn — a,b). The quantity Ym is the radiation pressure function given by 

4 ^ 

Yjn = —^^{n + l){a„i,n + arn,n+l +'2arn,naTn,n+l +'2l3rn,nl3rn,,n+l), m = a,b, (41) 

^™ n=0 

where Q!m,n and /3m, n are the real and the imaginary parts of Sm.n^ respectively. Moreover, when Atj = then 
Finally, the total radiation force exerted on the sphere by the dual-frequency plane wave is given by 

frad(0 = Trrg [EaYa + EbYb + ^A.;5R {^Ac.e-''^-*}) e,. (42) 



D. Numerical evaluation of the dynamic radiation force 

The dynamic radiation force function Yao; was evaluated numerically in Matlab 6.5 (MathWorks, Inc.). Two 
different materials were chosen in this evaluation: brass and stainless steel. The physical constants of the analyzed 
spheres are given in Table U The surrounding fluid of the sphere was assumed to be water with density po = 
1000 Kg/m'^, compressional velocity cq = 1500 m/s. The radius of the sphere is ro = 0.5 mm. 

We are interested here in analyzing how the dynamic radiation force changes with the center frequency /q. We 
evaluate the function Y/^^ as a function of /cq^o varying from 0.1 to 10. The difference frequency A/ was fixed to 
0, 50, 100 kHz. To assure that the center frequency /o is always positive, we assume that /o varies from 50 kHz to 
4.77 MHz. 

Before presenting the numerical evaluation results, let us take a closer look to the contribution of the second-order 
velocity potential to the dynamic radiation force given in Eq. (|29|l . For the frequency-range considered here, the 
energies Eq and E^^ are about the same order of magnitude. In this case, the numerical evaluation of i?o, see 
Eq. Ij40ll . for the chosen frequency-range has shown that the dimensionless amplitude of this quantity (10""^) is 
much smaller than the unit. In fact, for the frequency range used in the simulations we have Afcrg ^ 1. Therefore, 
we may neglect the contribution of i?o to the dynamic radiation force function given by Eq. ijifijl hereafter. 

In Figs. and 21 we see the magnitude of the dynamic radiation force function |Fal^|- The inset of the figures 
shows the phase of Yauj. We can see that when A/ — the function |1a;^ | of both materials is equal to the radiation 
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pressure function Ym, as expected. The dynamic radiation force function of both material exhibits a fluctuation 
pattern (dips and peaks) due to resonances of the ultrasound wave inside the sphere. The fluctuations depend 
on resonances of the scattering function Sm,n, which is related to the material parameters (density, compressional 
and shear speed of the wave) . No signiflcant changes in Yauj (magnitude and phase) is observed as the difference 
frequency varies from to 50 kHz. Further, the phase remains practically constant with zero value. This occurs 
because at 50kHz we have Afcro = 0.02, which implies kaVo ~ hro. Thus Yaui approaches to Ym- 

For the brass sphere, a prominent peak occurs in \Yaui\ with A/ = at fco'^o = 3.55, as seen in Fig. El When 
the difference frequency is 100 kHz, the peak change its position to 3.27 and the whole fluctuation pattern changes. 
However, the fluctuation form follows that one of A/ = with smaller amplitudes. The phase also presents 
fluctuation whose amplitude is approximately 7r/6rad. 

The function \Yauj\ for the stainless steel sphere with A/ — presents dips, as shown in Fig. 31 The flrst dips 
occurs at kgro = 5.17. At A/ = 100 kHz, the fluctuations in Yauj have a different pattern with smaller amplitude 
compared to the case of A/ — 0. The phase of Yauj is practically constant when A/ — 50 kHz, except when k^ro < 1. 
For A/ = 100 kHz, the phase exhibits fluctuations with amplitude of about 7r/6rad. The phase fluctuations follows 
the pattern exhibited in the magnitude of Yauj ■ 

Now we analyze how the dynamic radiation force changes as A/ varies. This is particularly useful to show possible 
resonances on the dynamic radiation force function of the spheres for a given center frequency /q. In Fig. the 
function IYai^I is plotted as A/ varies in from —100 kHz to 100 kHz. The phase of Ya^; is shown in the inset of the 
flgure. The frequency-range of A/ was chosen to reveal symmetry properties of Yalo- We flxed koro to 3.55 and 
4.82 for the dashed and solid lines, respectively. These values correspond to the flrst peak and dip, respectively, in 
Fig- El In both cases the IYac^I is symmetric. The phase of Yau is antisymmetric. Both magnitude and phase of 
Yaui change shape considerably as fegro changes. The concavities of \Yau,\ in Fig.Ofollow those shown in Fig.E|for 
the resonance points 3.55 (peak) and 4.82 (dip). 

The plot of the function Yaw of the stainless steel sphere as A/ varies from —100 kHz to 100 kHz is shown in 
Fig- El The inset of the flgure plots the phase of |1au|- The quantity fcoro was flxed at the flrst and second dips. 
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which corresponds to 5.17 and 6.85 for the dashed and soUd lines, respectively. The function IYai^I for both values 
of fcoro exhibited a symmetric concave pattern. The phase showed an antisymmetric pattern in both cases. The 
concavities of \Yauj\ seen in Fig. Elfollow those presented in Fig.^lfor points 5.17 and 6.85. 

The comparison of the normaHzed amplitude of the static and the dynamic radiation forces on the stainless sphere 
as the center frequency varies can be seen in Fig. [3 The radiation force functions Ya and Yt are given by Eq. Ij41|l . 
The difference frequency was fixed to 10 kHz. The amplitudes are normalized by the highest ultrasound energy 
density Ea times the cross-section area of the sphere. In this figure, the solid Hne stands for the amplitude of the 
static radiation force, Ya + (k^/kaYY},. The dotted line corresponds to {ki,/ka)\Y/:^^\. The magnitude of the both 
radiation forces are about the same. Though we changed the difference frequency in a range up to 100 kHz, the 
magnitude of the static and the dynamic radiation force remained approximately the same. 



IV. CONCLUSIONS 



We have presented a theory to calculate the dynamic ultrasound radiation force exerted on an object by a dual- 
frequency CW ultrasound beam in lossless fluids. The theory is valid for beams with any spatial distribution. The 
amplitude of the induced vibration by the dynamic radiation force on the object was assumed to be much smaller 
than its characteristic dimensions. No assumptions have been made on geometric shape of the object. A formula for 
the dynamic radiation force was obtained - see Eq. ll22ll - in terms of the first- and second-order velocity potentials. 
The dependence of the dynamic radiation force with the nonlinear parameter B/A of the medium was analyzed. 

We have calculated the dynamic radiation force exerted on a solid sphere by a dual-frequency CW plane wave in 
water. The dynamic radiation force is proportional to the cross-section area of the sphere, the dynamic ultrasound 
energy, and the dynamic radiation force function. The contribution of the first-order velocity potential to the 
radiation force, accounted by Eq. (jSOJ, was neglected because we considered that the dislocation of the sphere is 
very small. The contribution of the medium nonlinearity to the dynamic radiation force is negligible if A/cro ^ 1 
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FIG. 5: Dynamic radiation force function as A/ varies for the brass sphere. Legend: dashed line (fcoro = 3.55) and solid 
line (fcoro = 4.82). 

in a weak nonlinear medium. However, it may become more significant in strongly nonlinear mediums or when 
Afcro ^ 1. Numerical evaluation of the dynamic radiation force function revealed a fiuctuation pattern as the center 
frequency varies. The fluctuations are similar to those present in the static radiation force function. Analysis of the 
static and the dynamic radiation force on the sphere has shown that they have approximately the same magnitude. 
Experimental measurements of the magnitude of the static and the dynamic radiation force on a stainless steel 
sphere have conflrmed this prediction 

In conclusion, the presented dynamic radiation force formula ll22|l generaHzes Yosioka's formula which stands 
only for the static radiation force. The dynamic radiation force formula can be extended for a multi-frequency CW 
ultrasound beam. This is particularly useful to study pulsed radiation force in which the incident ultrasound beam 
can be expanded in time as a Fourier series. 
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APPENDIX A: SECOND-ORDER VELOCITY POTENTIAL 



To calculate the contribution of the second-order velocity potential to the dynamic radiation force on the sphere, 
we consider the lossless Burger's equation 
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FIG. 6: Dynamic radiation force function as A/ varies for the stainless steel sphere. Legend: dashed line (fcoro — 5.17) and 
solid line (fcoro = 6.85). 



where v is the velocity particle in the z-direction, e = 1 + B/{2A), and r = t — z/cq is the retarded time. The 
source wave form is given by 

v{0,t) = i;o[sin(cJar) + sin(ti;br)], 

where vq is the peak amplitude of the velocity particle at the wave source {z — 0). Hence, the second-order velocity 
particle at the difference frequency in the pre-shock wave range is given by 

2 

J2) _ _^_}!^i^kzBin{/\ut - Akz), 
2co 

This approximated solution is vaHd for evoAkz/co ^ 1. We obtain the second-order velocity potential at the 
difference frequency by integrating this equation over z. Thus, in complex notation, the amplitude of the velocity 
potential at Atu is 



Notice that the time-dependent integration constant was dropped because it evaluates zero in the closed surface 
integral of Eq. 



14 



1 r 




FIG. 7: Comparison on the magnitude of the static and the dynamic ultrasound radiation force on the stainless steel sphere 
with A/ = lOkHz. Legend: dotted line represents {kb/ka)\YAw\ and solid line stands for Ya + {kb/ka)^Yb. 
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